Document made available under the 
Patent Cooperation Treaty (PCT) 



International application number: PCT7US2005/0 10867 
International filing date: 30 March 2005 (30.03.2005) 

Document type: Certified copy of priority document 

Document details: Country/Office: US 

Number: 60/558,174 

Filing date: 31 March 2004 (31.03.2004) 



Date of receipt at the International Bureau: 21 May 2007 (21.05.2007) 



Remark: Priority document submitted or transmitted to the International Bureau in 
compliance with Rule 17.1(a) or (b) 




World Intellectual Property Organization (WIPO) - Geneva, Switzerland 
Organisation Mondiale de la Propriete Intellectuelle (OMPI) - Geneve, Suisse 



THIS IS TO CERTIFY THAT ANNEXED HERETO IS A TRUE COPY FROM 
THE RECORDS OF THE UNITED STATES PATENT AND TRADEMARK 
OFFICE OF THOSE PAPERS OF THE BELOW IDENTIFIED PATENT 
APPLICATION THAT MET THE REQUIREMENTS TO BE GRANTED A 
FILING DATE. 



APPLICATION NUMBER: 60/558,174 
FILING DATE: March 31, 2004 

RELATED PCT APPLICATION NUMBER: PCT/US05/10867 



THE COUNTRY CODE AND NUMBER OF YOUR PRIORITY 
APPLICATION, TO BE USED FOR FILING ABROAD UNDER THE PARIS 
CONVENTION, IS VS60I 558,17 4 



t :NI I I 1> STATES DKPARIMEY) OF COMMERCE 



United Sink's Patent ami Trademark Office 



May 15, 2007 




Under Secretary >d'< t«nimervc 
for Intellectual Property 
and Director of the United States 
Patent and Trademark Office 



Certified hy 




CO 
CO 



Please type a plus sign (+) inside this box ^ ["+] 

PTO/SB/16 (10-01) 
Approved for use through 10/31/2002. OMB 0651-0032 
Patent and Trademark Office; U.S. DEPARTMENT OF COMMERCE 
Under the Paperwork Reduction Act of 1 995, no persons are required to respond to a collection of information unless it displays a valid OMB control number. 

PROVISIONAL APPLICATION FOR PATENT COVER SHEET 

This is a request for filing a PROVISIONAL APPLICATION FOR PATENT under 37 CFR 1 .53 (c). 



Express Mail Label No. EL 623 309 525 US 



Fathi M. 
Khurram 



INVENTOR(S) 



Family Name or Surname 



Salem 
Waheed 



(City ai 



pr Foreign Country) 



Okemos, Michigan 
Lansing, Michigan 



□ Additional inventors are being named on the separately numbered sheets attached hereto 



TITLE OF THE INVENTION (280 characters max) 



MULTI-USER DETECTION IN CDMA SYSTEMS 



CORRESPONDENCE ADDRESS 



Direct all a 
[3 Customer Number 



Type Customer Number here 



Individual Name 



City 



Harness, Dickey & Pierce, P.L.C. 



Bloomfield Hills 



USA 



Telephone 248-641-1600 



Fax 248-641-0270 



ENCLOSED APPLICATION PARTS (check all that apply) 

EZ] 



j Specification Number of Pages 
3 Drawing(s) Number of Sheets | 4 [ 

3 Application Data Sheet. See 37 CFR 1.76 



□ CD(s), Number L 



3 Other (specify) | Postcard 



^ Specification Filed in English 



METHOD OF PAYMENT OF FILING FEES FOR THIS PROVISIONAL APPLICATION FOR PATENT (check one) 



£3 Applicant claims small entity status. See 37 CFR 1 .27. 
A check or money order is enclosed to cover the filing fees 

I | The Commissioner is hereby authorized to charge filing r— 
fees or credit any overpayment to Deposit Account Number: |_ 
□ Payment by credit card. Form PTO-2038 is attached. 



FILING FEE 
AMOUNT ($) 



The invention was made by an agency of the United States Government or under a contract with an agency of 
the United States Government. 

S No. 

□ Yes, the name of the U.S. Government agency and the Government contract number are: . 



Respectfully submitted, 
SIGNATURE 



TYPED or PRINTED NAME Jennifer S. Brooks 



JZ s/<L^^ Date i 3/31/04 i 

' QPrSICTOATIOM MD I cT 



TELEPHONE 248-641-1600 



REGISTRATION NO. 
(if appropriate) ' 
Docket Number: | 6550-000087 j 



USE ONLY FOR FILING A PROVISIONAL APPLICATION FOR PATENT 



Attorney Docket No. 6550-000087 
MULTI-USER DETECTION IN CDMA SYSTEMS 

FIELD OF THE INVENTION 
[0001] The present invention relates to multi-user detection in Code 
Division Multiple Access (CDMA) systems. 

BACKGROUND AND SUMMARY OF THE INVENTION 
[0002] Code Division Multiple Access (CDMA) is based on spread- 
spectrum technology and is a dominant air interface for the proposed modern 3G 
and 4G wireless networks. The transmitted CDMA signals propagate through noisy 
multipath fading communication channels before arriving at the receiver of the user 
equipment (UE). In contrast to classical single-user detection (SUD) computational 
techniques, which do not provide the requisite performance for modern high data 
rate applications, multi-user detection (MUD) approaches require a lot of a-priori 
information not available to the UE. 

[0003] The conventional detection schemes for CDMA signals only exploit 
second order statistics among user codes. Practically, however, the underlying 
user data symbol sequences are in general mutually (near-) "independent". This is 
a key assumption, which propels the application of info-theoretic learning 
approaches such as information maximization and minimum mutual information to 
the realm of CDMA. The use of these computational techniques is justified since a 
wide sense stationary slowly fading multipath CDMA environment can be 
conveniently represented as a linear multi-channel convolution model. The 
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received CDMA signal can be considered as a sum of several non-gaussian 
random variables generated by the linear convolutive transformations of statistically 
(near-) independent component user variables. This linear transformation accounts 
for the user spreading codes, the cell scrambling codes (in case of a cellular 
architecture), multiple channel paths and slowly fading channel effects. The 
present invention estimates a linear transformation to counteract, as "optimally" as 
possible, the effects of a first transformation -resulting in the recovery of the 
original user signals under the constraint of knowing only the user's signature code 
(and the corresponding cell scrambling code for a two stage implementation). 

[0004] Blind Source Recovery (BSR) is the process of estimating the 
original "independent" user-specific symbol sequences independent of, and even in 
the absence of, precise system/channel identification. In typical downlink signal 
processing, where many of the system parameters are unknown, including the 
codes for co-existing users at any instant of time, one can use the blind techniques 
for better estimate of the user-specific signals. Alternately, Blind Multi-User 
Detection (BMUD) computational techniques based on the Natural Gradient Blind 
Source Recovery (BSR) techniques in both feedback and feedforward structures 
can be used. The "quasi-orthogonality" of the spreading codes and the inherent 
"independence" among the various transmitted user symbol sequences form the 
basis of the proposed BMUD methods. The inventive structures and computational 
techniques demonstrate promising results as compared to the conventional 
techniques comprising, e.g., Match Filter (MF), RAKE and block LMMSE. The 
inventive computational techniques can be implemented either using the batch or 
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the more computationally efficient instantaneous update methods. Although batch 
implementations exhibit better performance, it is however accompanied by longer 
latency and more involved implementation structures not suitable for a UE/MS. 
Focus is placed on the instantaneous (or on-line) performance of the BSR 
computational techniques only, which still exceeds the performance of other 
approaches. 

[0005] This (on-line) detection technique can be easily extended to other 
CDMA implementations, using relatively short scrambling codes, but becomes 
impractical in WCDMA downlink using long scrambling codes. In spite of the fact 
that very low bit error rates (BER) can be achieved with this technique and the 
detection process does not even require the knowledge of user's own signature 
code, the recovered signal stream is at the symbol level with no user identification. 
Further, inherent sign and permutation ambiguities exist in BSR (scaling is not 
relevant as the recovered streams are either BPSK or QPSK etc.); therefore user 
identification is not possible unless some preamble or pilot data is transmitted 
periodically. This periodic requirement stems from the dynamic nature of the 
wireless communication scenario where users may dynamically enter or exit the 
system. The environment structure also varies widely due to the mobility of the 
MS/UE. 

[0006] With these practical constraints in mind, new computational 
techniques are proposed by an infusion of info-theoretic learning computational 
techniques such as static Blind Source Recovery (BSR) (or Independent 
Component Analysis, ICA) and Principal Component Analysis (PCA) in the existing 
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structure of a RAKE receiver. The purpose of this additional info-theoretic stage is 
to counter, as best as possible, the unmodeled multiple access interference (MAI) 
and the additive (white Gaussian) noise in the channel. Further, use of a simple 
info-theoretic stage does not make the receiver structure too complex (in fact, it is 
simpler than most other proposed adaptive LMMSE implementations. RAKE-PCA 
uses up to second order statistics, as compared to RAKE-BSR, which utilizes 
higher order statistics. This results in slightly simpler update structure for the 
RAKE-PCA, but the performance of the RAKE-BSR is found to be better than 
RAKE-PCA. Further, assuming the score-function for the ICA update law to be 
chosen properly, the resulting equalization matrix in case of RAKE-BSR has 
relatively smaller element values (energy) as compared to the corresponding matrix 
for RAKE-PCA, which can be translated to the need of fewer memory bits. Lastly, 
both RAKE-BSR and RAKE-PCA use all the available user information, so 

that there are no issues of user identification and sign ambiguities in this case. The 
main advantage of both the RAKE-BSR and RAKE-PCA computational techniques 
is the improved BER performance for the UE/MS without the need of any additional 
information than what a standard RAKE receiver already has. The proposed 
computational techniques can be applied to both generic DS-CDMA and modern 
3G UMTS implementations. 

[0007] Further areas of applicability of the present invention will become 
apparent from the detailed description provided hereinafter. It should be 
understood that the detailed description and specific examples, while indicating the 
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preferred embodiment of the invention, are intended for purposes of illustration only 
and are not intended to limit the scope of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0008] The present invention will become more fully understood from the 
detailed description and the accompanying drawings, wherein: 

[0009] Figure 1 is a block diagram illustrating a signal generation model for 
a QPSK DS-CDMA system; 

[0010] Figure 2 is a block diagram illustrating a feedforward demixing 
structure in accordance with a first embodiment of the present invention; 

[0011] Figure 3 is a block diagram illustrating a feedback demixing 
structure in accordance with a second embodiment according to the present 
invention; and 

[0012] Figure 4 is a is a block diagram illustrating a feedback demixing 
structure in accordance with a third embodiment according to the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0013] The following description of the preferred embodiments is merely 
exemplary in nature and is in no way intended to limit the invention, its application, 
or uses. 

[0014] The present invention includes three embodiments of demixing 
structures providing new MUD detection systems and methods, and two new types 
detectors using BSR techniques. The new MUD detection systems and methods 
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are discussed as Natural Gradient Blind Multi-User Detection (BMUD) 
computational systems and methods. The two new detectors are RAKE-Blind 
Source Recovery (RAKE-BSR) and RAKE-Principal Component Analysis (RAKE- 
PCA) Detectors. These detection systems, methods, and detectors are discussed 
with reference to a convenient convolutive signal model representation of DS- 
CDMA systems discussed with reference to Figure 1. 

[0015] In a typical downlink synchronous DS-CDMA system employed for 
indoor ATM and certain ad-hoc wireless networks, each user's data 10 is spread 
using an individual signature waveform (or spreading code), then the data 10 for all 
users is combined and transmitted over multipath AWGN channel 12 by the Base 
Station (BS) 14. Each User Equipment (UE) or Mobile station (MS) synchronizes 
itself with the BS using the broadcast synchronization/pilot channels, once 
synchronized the BS and UE/MS can communicate on the traffic channel 
(comprising of both data and control streams). Assuming the data transmission to 
be QPSK, i.e., comprising of two composite data channels created by a serial-to- 
parallel (S/P) stage, which are constellated in quadrature. At the UE/MS receiver, 
the received signal is first passed through a chip-matched filter and sampled at chip 
rate. 

[0016] Considering a total of K active users in an L multipath environment 
and N transmitted symbols during the observation frame T F , the received signal is 
given by 

KO-f f S V^*.(«K('K('-»T-r,) + n(,) <1) 
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where ^is the energy of the n' h symbol for the ^user, fc t (n)e{±l±»}is the n' h 
complex symbol for the it" 1 user, fyand r, are the /""path's gain co-efficients and 
delay, respectively. n(t) is the additive noise and s k (t) is the k' h user's signature 
code (or spreading sequence) generated by 

G-l 

s k (t) = Y J C( k {m)p(t - mT c ); a k (m)e {- 1,1}, 0 < m < G - 1 (2) 

m=0 

a k {m) is a real spreading sequence (i.e., any of the standard CDMA PN codes, 
such as the Gold, Walsh-Hadamard, Kasami etc.) sequence for the k' h user 
containing G chips per symbol, i.e., G = T 6 /T C .p(t) is a chipping pulse of duration 
T c .T fc being the symbol period. 

[0017] Under the assumption of time-invariance, the model (1) can be 
more compactly written in a vector-matrix format at 

r=HSb + n (3) 
where, H is a (NG + L-l)xNG multipath propagation co-efficient matrix containing 
the channel coefficients. S is a NG xNG block diagonal matrix with the matrix of 
spreading codes forming the diagonal elements, b is an NG -d vector containing 
the user symbols, while n is the (NG + L-l)-d channel noise vector with 
covariance matrix Q. The structure of the above defined matrices and vectors is 
given by 
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|i ••. ; | 
'•• o | 
|o ■•• K I 
f ••• ••• •■ I 

L° ••• °vJ 

S=diag[sS --s\,S = [s, s 2 ■•■ sj 

b = [b{l) T b(2) T ■ ■ ■ b{N) T \, b(n) = [•/ 7 I7l(n) V £ 2n 6 l(n) ■ ■ ■ Je^b K {n)\ 
[0018] The compact linear model (3) is useful in deriving the linear 
detectors such as matched filter (MF), linear minimum mean squared error 
(LMMSE) etc. for recovery of the transmitted symbol train for a desired user. 
However, the primary disadvantage of this model is the prohibitive dimensions of 
the constituent matrices, especially with longer frame durations and larger 6, the 
matrices become excessively large, making this model unsuitable for any real-time 
implementation at UE/MS. 

[0019] Alternately, the signal model can be represented as a linear 
convolutive model, i.e., during symbol time, the received chip data is constituted of 
the chips corresponding to the currently transmitted symbol and its delayed 
multipath components as well as delayed chips from some previously transmitted 
symbols and the channel added noise. In this formulation, G chips arriving at the 
UE/MS during the n' h symbol time are computed as the sum of the chips from L 
multipaths of the n' h transmitted symbol and the multipath components of the 
previous J- 1 symbols (n-l,...,n- J -l), where 
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(4) 



t l being the maximum chip delay in L multipaths. The n' h received symbol can be 
expressed as 

r n (t) 4>k MfetifZ h ('k (< - "T- T, )+ n n 

Z 2>* ( n - A/* (" - I>< ('K (' - (» - ^ - 7, ); "T < r < (n + 1)T (5) 

j=l k=l 1=0 

[0020] Under the assumption L<G, the above model can be 
expressed just in terms of the current and the preceding symbol. That is, the 
multipaths with delay greater than symbol period either do not exist or are weak 
enough to be ignored. In this case, the output samples of the chip-matched filter 
can be written as: 



k=\ L '=0 /=0 J 

where, z w and z w are G-d early and late code vectors, i.e., 



s k [G-r,+l} ••• sJg]0 --0 



(6) 



(7) 



(8) 



t, is the discretized delay satisfying the constraint 0<r, <T 6 . Imposing time 
invariant constraints, the multipath slowly fading environment model (6) can be 
represented in the form 

r n =li 0 b n +U x b n _ l +n n (9) 

9 
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where b n and b n _ x are the K-d vectors of current and previous symbol for all the 
K users. H 0 and H,are GxK mixing matrices with the structure 

H 0=L H 0,0 H 0.1 H 0,K J> H 1 =L H 1,0 H l.l H 1.kJ 

such that 

#o,*=V^I>/Z« (10) 

/=0 

H, t =v^i;^z a (id 

[0021] and e 0 e x represent the energy of the current and the previous 
symbol respectively. 

[0022] The Natural Gradient Blind Multi-User Detection (BMUD) system 
and method is discussed below with reference to Figures 2-4. Blind Multi-user 
Detection (BMUD) is the process to blindly estimate all the user symbol sequences 
directly from the received composite CDMA signal using the Blind Source Recovery 
(BSR) techniques. BSR framework implies recovery of original signals from 
environments that may include transient, convolution and even non-linearity. The 
linear BSR computational techniques have been developed for linear convolutive 
mixing environments by the minimization of mutual information (Kullback Lieblar 
Divergence) using the natural gradient subject to the structural constraints of a 
recovery network. The natural gradient BMUD network can be either in the 
feedforward or the feedback configuration, in which the proposed BMUD system 
and method adaptively estimates a set of matrices to counter the linear convolutive 
environment model (9). In each case, an input receiveseither the linear convolutive 
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environment model r n or a whitened version r n w of the linear convolutive 
environment model r n . Parametric matrices 16A-16C W 0 and W k are adaptable to 
estimate independent user symbols y n at an n' h instant based on the linear 
convolutive environment model r n or the whitened version r™ of the linear 
convolutive environment model r„ . A decision stage 18 interprets y n and 
estimates corresponding user symbol estimates b n also at the n' h instant. 

[0023] The justification for BMUD computational techniques is based 
on the convenient convolutive signal model representation of DS-CDMA systems, 
see (9), and the reasonable assumption that the various transmitted user symbol 
sequences are mutually "independent" as they are generated independently of each 
other. In this framework, both the transmitted sequence and the mixing matrices in 
the model (9) are unknown to the user. The only known entity to the user is the 
self-identification code. Other available prior information is the nature of transmitted 
data, which is assumed to be quaternary sub-guassian distribution, e.g., QPSK data 
distorted by the mutipath channel and additive noise. There exists enough 
information to apply the info-theoretic natural gradient Blind Source Recovery 
(BSR) system and method for BMUD in this case. 

[0024] Further assume that the DS-CDMA channel is not over-saturated 
and K < G . The proposed BMUD computational techniques do not require any pre- 
whitening of received data. However, in most modern WCDMA, G is chosen to be 
very large and in general K < G . Therefore, it is computationally advantageous to 
pre-process received data for dimension reduction to K which is the actual number 
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of principal independent symbol sequences (or users) in the received data. The 
process of pre-whitening will also remove the second order dependence among the 
received data samples and some of the additive noise. The data pre-whitening can 
be achieved either online using adaptive principal component analysis (PCA) 
computational techniques or it may be done Using an algebraic PCA estimate over 
a large batch (say N samples) of received data, i.e., 

R = [r t r 2 ■■• r N-l r N] 

with the correlation matrix 

A r =-l—RR T 
c N-l 

Then the whitening is achieved using the filtering matrix 

W = D~ U2 V T 
where 

D: represents the K-dim matrix of principle eigenvalues of the data correlation matrix 

Ac- 

V: represents the KxM matrix of principle eigen vectors of the data correlation matrix 

A c - 

[0025] The whitened version of (9) is given by 

ijT = W(H 0 b n + + »») = H 0 b n + Hfa-i (12) 

where r™ represents the G-Jwhitened data received during the n ,h symbol time 
and H 0 ,H, represent the equivalent GxK convolutive mixing matrices for the 
current and the previous symbol. 
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[0026] BMUD computational techniques blindly adapt a set of matrices to 
estimate the independent user symbols y n at the n' h instant. This is followed by a 
decision stage to interpret, as best as possible, y n and estimate the corresponding 
user symbol estimates b n also at the n' h instant. 

k=¥iy n ) (13) 
where y/{) \ represents the (nonlinear) decision stage. 

[0027] For the feedforward configuration discussed with reference to 
Figure 2, the BMUD stage output is computed as 

y n =wvZ+hvkC k 04) 
k=\ 

[0028] The update laws for this feedforward structure can be derived. For 
the feedforward parametric matrices W 0 and W k , the update laws are 

AW 0 ~{l-<p(y n )yZ )w 0 (15) 

and 

W k ~[l-(p{y n )y»y k -<p(y n )(r™ k f (16) 
where <p{) is an element-wise acting score function, and lis a K-d identity 
matrix. 

[0029] For the initialization of the computational technique, W 0 is chosen 
to be either an identity or a diagonally dominantly matrix, while all other matrices 
W k are initialized to have either random elements with a very small variance or as a 
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matrix of all zeros. Note that no matrix inversion is required for the feedforward 
computational technique. 

[0030] In a second embodiment discussed with reference to Figure 3 and 
hereinafter referred to as feedback configuration I, the output is estimated by 



V k=l J 

The update law for this structure using the natural gradient are given for the matrix 

W 0 by 



The matrices in this case are also initialized in a fashion similar to the feedforward 
case. However, note that at least one matrix inversion is required in this 
formulation. 

[0031] In a third embodiment discussed with reference to Figure 4 and 
hereinafter referred to as feedback configuration II, the feedback configuration 
implements the feedback structure without the need for any matrix inversion. For 
this feedback configuration II, the BMUD stage output is computed as 





(18) 



While for the feedback matrices W k , the update law is 




(19) 



K 



yn=w 0 rf-Zw k y n - k 



(20) 



k=l 



The update laws for this feedback structure are given by 




(21) 



14 



Attorney Docket No. 6550-000087 

&W k <*\[-rty a )y*)tiy k +ay n )yn_ k (22) 
[0032] In case the channel can be estimated, the performance of the 
proposed BMUD computational techniques may be adjudged by the diagonalization 
of the absolute value of the global transfer function. The global transfer function 
presents the combined effect of the complex mixing and demixing transfer 
functions. For the two symbol convolutive models for the case of L<G, the global 
transfer function for the natural gradient computational techniques in the z-domain 
are given by: 

G = G 0 +G lZ ~ l +G 2 z~ 2 (23) 
where, for the feedforward configuration 
G 0 =W 0 H 0 =W 0 WH 0 , 

G l = Wb#i + WiH 0 = W 0 WH l + W X WH 0 and (24) 

G 2 =W 1 H l =W l WH l 

while, for the feedback configuration I 

G 0 =Wo 1 H 0 =Wq 1 WH 0 , 

Gi = Wq 1 (H l -W l ) = Wq 1 {WH x - Wi ) and (25) 
G 2 =0 

and for the feedback configuration II 
G 0 =W 0 H 0 =W 0 WH 0 , 

Gi = WqHy +Wi= W 0 WH l + Wi and (26) 
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G 2 =0 



[0033] The proposed BMUD computational techniques as formulated 
result in recovery of the user symbols directly. The computational techniques can 
be conveniently applied to DS-CDMA systems using only user-specific spreading 
sequences. They may also be extended to CDMA systems using relatively short 
scrambling codes, though the dimension o matrices may still become large for DSP 
implementations in a UE/MS. The WCDMA system uses long codes in the 
downlink, making the application of these computational techniques impractical 
because of the requisite dimension of the demixing network matrices. 

[0034] The RAKE-Blind Source Recovery (RAKE-BSR) and RAKE- 
Principal Component Analysis (RAKE-PCA) Detectors are now described. RAKE- 
BSR and RAKE-PCA are two new proposed adaptive detectors, which utilize the 
same knowledge as a RAKE receiver. An info-theoretic adaptive weighting matrix of 
dimension GxG is introduced into the RAKE structure, which gives a big 
performance boost to the RAKE receiver. The performance of RAKE-BSR/RAKE- 
PCA exceeds the performance of LMMSE detectors under the conditions of high 
network congestion, imprecise channel estimation, and unmodeled inter-cellular 
interference etc.. The closed form structure of these proposed detectors is given by 



where W = diag[AA-A], and A is the GxG matrix that is adaptively estimated 
either using static BSR (ICA) or PCA computational techniques. 



b. 




(27) 
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[0035] It is proposed to adapt the matrix A using the natural gradient 
update laws. However, there exist several other methds for ICA/PCA and any other 
suitable method may be used for these adaptations. This blind adaptation of the A 
matrix has several advantages and improves the performance of the overall 
equalization process in several ways. Firstly, it can counter artifacts in the 
estimated channel co-efficients H . Secondly, the channel estimation process (as in 
RAKE receivers) may be limited by the structure (such as number of fingers) and 
may estimate only a few of the dominant channel parameters. W stage tends to 
counteract this anomaly, as best as possible, and provides better performance than 
LMMSE in such cases. Thirdly, this adaptive stage minimizes the effect of the 
additive channel noise, which may have an unmodeled intricate underlying 
structure. Lastly, the natural gradient ICA/PCA computational techniques inherently 
reduce near-far problems by removing any ill conditioning in the signal space for all 
the users in the system. This results in all the mobile users in the system to have a 
BER performance similar to the average BER performance of the downlink channel. 
The matrix A is adaptively estimated using the update laws 



and <p{) is a nonlinear score function which depends on the underlying distribution 
structure of the signals involved. For QPSK signal, a suitable score-function is 



A(k + l) = A(k) + rj k AA(k) 



(28) 



where 




for static BSR (or ICA) 
forPCA 



(29) 



<Pi(Ji) = - «;(tanh(# Re{y f })+ tanh(# Im{y J)) 



(30) 
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[0036] Of these proposed RAKE-BSR/RAKE-PCA structures, RAKE-BSR 
exhibits relatively faster and more stable convergence. However, in CDMA 
systems the underlying code structure is chosen to be "orthogonal", and RAKE- 
PCA may converge to a slightly lower BER solution if the channel impairments are 
linear in nature. Note that in (27), if the channel estimate H is either not available 
or changes very dynamically, the detector can be estimated without using the 
channel estimate and the structure reduces to Matched Filter BSR/PCA, i.e., 



The performance of this structure is better than Matched Filter alone, and 
approaches RAKE performance as the underlying matrix A converges. 

[0037] The description of the invention is merely exemplary in nature and, 
thus, variations that do not depart from the gist of the invention are intended to be 
within the scope of the invention. Such variations are not to be regarded as a 
departure from the spirit and scope of the invention. 




for DS-CDMA Systems 
forWCDMA Systems 



(31) 



18 



Attorney Docket No. 6550-000087 

CLAIMS 

What is claimed is: 

1. A natural gradient Blind Multi User Detection (BMUD) network system 
that adaptively estimates a set of matrices to counter a linear convolutive 
environment model r n , the system comprising: 

an input receptive of at least one of the linear convolutive environment 
model r n and a whitened version r™ of the linear convolutive environment model r n ; 

parametric matrices W 0 and W k adaptable to estimate independent 
user symbols y n at an n' h instant based on at least one of the linear convolutive 
environment model r n and the whitened version r n w of the linear convolutive 
environment model r n ; and 

a decision stage interpreting y n and estimating corresponding user 
symbol estimates b n also at the n' h instant. 

2. The system of claim 1, wherein the system is networked in a 
feedforward configuration. 

3. The system of claim 2, wherein the decision stage is adapted to 
compute y n according to: 

y n =Wo<+jy k c k . 



19 



Attorney Docket No. 6550-000087 



4. The system of claim 2, wherein the parametric matrices W 0 and W k 
have update laws according to: 



where <p(-) is an element-wise acting score function, and I is a K-d identity 
matrix. 

5. The system of claim 2, wherein W 0 is initially chosen to be at least one 
of an identity and a diagonally dominantly matrix, while all other matrices W k are 
initialized to have at least one of random elements with a very small variance and as 
matrices of all zeros. 

6. The system of claim 1 , wherein the system is networked in a feedback 
configuration. 

7. The system of claim 6, wherein the decision stage is adapted to 
compute y n according to: 
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8. The system of claim 6, wherein the parametric matrices W 0 and W k 
have update laws according to: 

AW 0 oc-W 0 [l-<p{y n )y%); and 
AW k ocW 0 (q>(y n )yff_ k ), 
where <p{) is an element-wise acting score function, and I is a K-d identity 
matrix. 

9. The system of claim 1 , wherein the system is networked in a feedback 
configuration without need for any matrix inversion. 

10. The system of claim 9, wherein the decision stage is adapted to 
compute y n according to: 

yn=^r:-JW k y n _ k . 
k=l 

11. The system of claim 9, wherein the parametric matrices W 0 and W k 
have update laws according to: 

AW 0 oc(/-^ n )^)^ 0 ;and 

AW k - (/ - <p{y n ) y n )w k + <p(y n )y%_ k , 
where q>{) is an element-wise acting score function, and I is a K-d identity 
matrix. 
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12. The system of claim 1, further comprising a whitening filter 
preprocessing received data for dimension reduction to K, which is an actual 
number of principal independent symbol sequences in the received data, and to 
remove second order dependence among received data samples and additive 
noise. 

13. The system of claim 12, wherein the whitening filter whitens data 
online using adaptive principle component analysis computational techniques. 

14. The system of claim 13, wherein the whitening filter whitens data using 
an algebraic PCA estimate over a large batch of received data including N samples 
according to: 

R = [r { r 2 ••• r N-l r Nj 

with a data correlation matrix 

A r =— !— RR T . 
c N-1 

15. The system of claim 14, wherein the filter achieves the whitening using 
a filtering matrix according to: 

W = D- 1,2 V T 

where D represents a K-dim matrix of principle eigenvalues of the data correlation 
matrix A c , and V represents a KxM matrix of principle eigen vectors of the data 
correlation matrix A c . 
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1 6. The system of claim 1 2, wherein the filter is adapted to calculate the 
whitened version r" of the linear convolutive environment model r„ according to: 

r n w = W(H 0 b n + + n n ) = H 0 b n + tf^l - 

and the linear convolutive environment model r n is represented according to: 

where b n and b n _ x are the K-d vectors of current and previous symbol for all the 
K users, H 0 and H,are GxK mixing matrices with the structure 

H 0 =[Ho.o H oj - H m - H,,J SUCh that H^ k =V^2^^ • 

/=o 

H /,t = yR\Y, h ilki > and £ o, £ \ represent the energy of the current and the previous 

(=0 

symbol respectively. 

17. An adaptive detector utilizing knowledge utilized by a RAKE receiver, 
comprising: 

an adaptive weighting matrix introduced into a RAKE structure, 
wherein the matrix is adaptively estimated using at least one of Principal 
Component Analysis (PCA) computational techniques and static Blind Source 
Recovery (BSR) computational techniques based on Independent Component 
Analysis (ICA). 
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18. The detector of claim 17, further comprising a closed form RAKE 
structure according to: 



where W = diag[A A ••■ A ], and A is the matrix. 

19. The detector of claim 18, wherein the matrix is adapted according to 
natural gradient update laws. 

20. The detector of claim 19, wherein the matrix is adaptively estimated 
using natural gradient update laws according to: 

A(k + l) = A(k) + rj k M(k) 



and <p{) is a nonlinear score function which depends on an underlying distribution 
structure of involved signals. 

21 . The detector of claim 20, wherein the score function is: 




for DS-CDMA Systems 
forWCDMA Systems 



where 




for static BSR (or ICA) 
forPCA 



<Pi(yi) = Vi 



>i y t - a t (tanh(# Refo }) + tanh(# Imfr t })) . 
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22. The detector of claim 17, wherein a channel estimate H is at least one 
of not available and changes dynamically, and the detector is estimated without 
using the channel estimate, such that the structure reduces to Matched Filter (MF) 
BSR/PCA according to: 

_\spWr for DS-CDMA Systems 

i,MF-BSRIPCA ~ | S p C Hwr forWCDMA Systems 

23. A natural gradient Blind Multi User Detection (BMUD) method that 
adaptively estimates a set of matrices to counter a linear convolutive environment 
model r n , comprising: 

receiving at least one of the linear convolutive environment model 
r n and a whitened version r™ of the linear convolutive environment model r n ; 

adapting parametric matrices W 0 and W k to estimate independent user 
symbols y n at an n' h instant based on at least one of the linear convolutive 
environment model r n and the whitened version r„ w of the linear convolutive 
environment model r n ; and 

interpreting y n and estimating corresponding user symbol estimates 
b n also at the n' h instant. 

24. The method of claim 23, further comprising employing a feedforward 
network configuration. 
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25. The method of claim 24, further comprising computing y n according to: 



26. The method of claim 24, further comprising updating the parametric 
matrices W 0 and W k via update laws according to: 



where is an element-wise acting score function, and I is a K-d identity 
matrix. 

27. The method of claim 24, further comprising: 

initializing W 0 to be at least one of an identity and a diagonally 
dominantly matrix; and 

initializing all other matrices W t to have at least one of random 
elements with a very small variance and as matrices of all zeros. 

28. The method of claim 23, further comprising employing a feedback 
network configuration. 
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29. The method of claim 28, further comprising computing y n according to: 



30. The method of claim 28, updating the parametric matrices W 0 and W k 
via update laws according to: 



where <p{) is an element-wise acting score function, and I is a K-d identity 
matrix. 

31. The method of claim 23, further comprising employing a feedback 
network configuration without need for any matrix inversion. 

32. The method of claim 31 , further comprising computing y n according to: 






K 
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33. The method of claim 31, further comprising updating the parametric 
matrices W 0 and W k via update laws according to: 



where <p(-) is an element-wise acting score function, and I is a K-d identity 
matrix. 

34. The method of claim 23, further comprising preprocessing received 
data for dimension reduction to K, which is an actual number of principal 
independent symbol sequences in the received data, and to remove second order 
dependence among received data samples and additive noise. 

35. The method of claim 34, further comprising whitening data online using 
adaptive principle component analysis computational techniques. 

36. The method of claim 35, further comprising whitening data using an 
algebraic PCA estimate over a large batch of received data including N samples 
according to: 





R = r. r 2 



[r,r 2 --- r N-l r Nj 



with a data correlation matrix 




RR r . 
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37. The method of claim 35, further comprising employing a filtering matrix 
according to: 

W = D- U2 V T 

where D represents a K-dim matrix of principle eigenvalues of the data correlation 
matrix A c , and V represents a KxM matrix of principle eigen vectors of the data 
correlation matrix A c . 

38. The method of claim 34, further comprising calculating the whitened 
version r™ of the linear convolutive environment model r n according to: 

r B w = W(H 0 b n + Hfa-i +n n )^ H 0 b n + tf^-l - 
wherein the linear convolutive environment model r n is represented according to: 
r n =H 0 &„ +H,fc„_ 1 +n n 

where b n and b n _ x are the K-d vectors of current and previous symbol for all the 
K users, H 0 and H,are GxK mixing matrices with the structure 

H 0 =[H 0 ,o H w - Ho.J.H^Lh,,,, H,,, ... H U J such that H w =V^ZM« > 

(=0 

H lk =4e^Y J h l z kl , and s Q e x represent the energy of the current and the previous 
symbol respectively. 
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39. An adaptive detection method, comprising: 

introducing an adaptive weighting matrix into a RAKE structure, 
wherein the matrix is adaptively estimated using at least one of Principal 
Component Analysis (PCA) computational techniques and static Blind Source 
Recovery (BSR) computational techniques based on Independent Component 
Analysis (ICA). 

40. The method of claim 39, further comprising employing a closed form 
RAKE structure according to: 



where W = diag[A A • • • A ] , and A is the matrix. 

41. The method of claim 40, further comprising adapting the matrix 
according to natural gradient update laws. 
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42. The method of claim 41, further comprising adaptively estimating the 
matrix using natural gradient update laws according to: 
A(k + l) = A(k) + rj k M(k) 

where 



AA(k) = 



(I-y(k)y(k) H )A(k) forPCA 



and is a nonlinear score function which depends on an underlying distribution 
structure of involved signals. 

43. The method of claim 42, further comprising employing a score function 
according to: 

<Pi 0>i ) = v i yt ~ a i (tanh(A- Refo }) + tanh(# Im{y f })) . 



44. The method of claim 39, wherein a channel estimate H is at least one 
of not available and changes dynamically, the method further comprising estimating 
a detector without using the channel estimate, such that the detector structure 
reduces to Matched Filter BSR/PCA according to: 

hpWr for DS-CDMA Systems 

iMF-BSRIPCA~y S Hc H Wf for WCDMA Systems 
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ABSTRACT OF THE DISCLOSURE 
A natural gradient Blind Multi User Detection (BMUD) network system and 
method adaptively estimates a set of matrices to counter a linear convolutive 
environment model. Feedforward and feedback network structures may be 
implemented, with or without matrix inversion. In other aspects, an adaptive 
weighting matrix is introduced into a RAKE structure, and the matrix is adaptively 
estimated using Principal Component Analysis (PCA) computational techniques 
and/or static Blind Source Recovery (BSR) computational techniques based on 
Independent Component Analysis (ICA). 
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